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Abstract. We compute the 0(1/A^^) correction to the critical exponent 2A 
= — f3'{gc) for the chiral Gross Neveu model in arbitrary dimensions by sub- 
stituting the corrections to the asymptotic scaling forms of the propagators 
into the Schwinger Dyson equations and solving the resulting consistency 
equations. 
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1 Introduction. 



The large N expansion has proved to be a useful technique in providing an 
alternative method of analysing renormalizable quantum field theories which 
possess N fundamental fields. For example, it was demonstrated in jl| that 
the (perturbatively massless) two dimensional Gross Neveu models with dis- 
crete and continuous chiral symmetry exhibit dynamical mass generation as 
well as being asymptotically free field theories with this technique. This for- 
mer property cannot be accessed in conventional perturbation theory and 
such toy models have proved to be important as a step to understanding 
similar phenomena in complicated four dimensional gauge theories such as 
quantum chromodynamics. One of the limitations, however, of a conven- 
tional large analysis is that one cannot proceed much beyond the leading 
order due in part to the appearance of intractable integrals whose treatment 
is hindered primarily by the presence of the dynamically generated mass. 

Methods to overcome this problem have been developed in ||2|, ^ and ap- 
plied initially to the 0{N) bosonic a model to solve that theory to 0{1/N'^). 
The technique exploits the properties of the field theory at its d-dimensional 
critical point which is defined to be the non-trivial zero of the /3-function. 
The simplifying feature of studying a model there is the additional symme- 
try present. As the /3-function vanishes the Green's function of the fields 
are conformally invariant and obey simple power law behaviour which, of 
course, has a massless form. In essence one is using the conformal bootstrap 
approach to solve the field theory, Q. As a consequence, one avoids one of 
the difficulties of the conventional large N approach in that the previously 
intractable integrals become calculable using techniques developed to com- 
pute massless Feynman diagrams, ||5|. This technique, known as uniqueness, 
is in effect a method of conformal integration in d-dimensions. Essentially 
when the sum of the exponents of bosonic propagators at a 3- vertex in co- 
ordinate space sum to the dimension of space time, the vertex is replaced 
by a triangle of propagators whose exponents are related to the original ver- 
tex. This represents an integration within the overall Feynman diagram. In 
essence it is the Yang Baxter star triangle relation, but is more general since 
it is valid in arbitrary (fixed) space time dimension. Moreover, as one is now 
dealing with a scale invariant situation the problem of solving the model be- 
comes one of computing the critical indices or exponents of the fields which 
will have an anomalous portion in addition to the canonical dimension. The 
anomalous piece is the part which carries the quantum properties of the 
model and from the universality criterion of statistical mechanics (see, for 
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example, Q), it will be a function of the spacetime dimension and the basic 
parameters of the model, such as N in the case of a large N expandable 
theory, and can be calculated order by order in The equations which 

one uses to do this are deduced by representing the Dyson equations of the 
Green's functions in the critical region of the theory from which one can 
deduce self consistent equations whose solution yields the exponents, [§]. 
The analytic expressions which result for the anomalous dimensions can be 
related to the respective renormalization group functions precisely at the 
critical point, [§, |^. One important function which characterizes a field 
theory is the /^-function which governs how the coupling constant behaves 
with the renormalization scale. Ordinarily one calculates it order by order 
in powers of the perturbative coupling constant, which is assumed to be 
small. In the large critical point approach which we concentrate on here 
it can, however, be determined by computing the critical exponent 2A = — 
P'idc) where gc is the critical point. Whilst this may seem to be an indirect 
way to proceed it is important to note that the 0{1/N'^) corrections to the 
/3-function of the 0{N) bosonic a model could only be deduced in this way, 

i- 

Following the pioneering work of ||2|, 0] the method has been extended to 
analyse models with fermions, [7-11]. In particular the 0{N) Gross Neveu 
model has been solved at O^l/N"^) with the fermion and auxiliary field 
anomalous dimensions, [0, P], and the /3-function exponent A, fo], 10|, all 



having been computed in arbitrary dimensions. Now that matter fields have 
been successfully incorporated within the critical point approach it is possi- 
ble to examine other models to the same precision. One such model which is 
currently relevant and related to the 0{N) Gross Neveu model is the chiral 
Gross Neveu model, [|l|, |l2|, which possesses a continuous chiral symmetry. 
This was first discussed in by Nambu and Jona-Lasinio in four dimen- 



sions and later in [13| where the connection with hadronic physics was first 
introduced. Moreover, four fermi interactions have recently been a subject 
of intense study as it provides a possible alternative to the Higgs mechanism 
in the standard model. There a fermion fermion bound state plays the role 
of the Higgs boson, |14, ^] and so it is important to have as complete a 



picture as possible of the quantum structure of the model. Therefore, the 
aim of this paper is to compute the exponent A of the chiral Gross Neveu 
or Nambu-Jona-Lasinio model at 0{1/N'^) which relates to the /3-function. 
This is one order beyond what has been computed before in this model. 
One property that the chiral Gross Neveu model shares with other models 



in two dimensions is that it possesses an exact S'-matrix, |16], from which 
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an expression for the exact mass gap has recently been deduced, |17|. The 
S-matrix has been constructed on the assumption that 3-particle S'-matrix 
elements factorize into the product of the constituent 2-particle elements 
which is a consequence of the underlying Yang Baxter triangle relation as- 
sociated with integrable field theories. Whilst the model is asymptotically 
free and therefore possesses the property of dimensional transmutation in 
two dimensions, it might appear that one has already determined the (3- 
function of the 2-dimensional model. However, such a statement overlooks 
the implicit nature of the situation. As far as we are aware, the explicit form 
of the 2-dimensional /3-function which can in principle be deduced from the 
S'-matrix spectrum, has never been achieved in detail. In providing new in- 
formation on the 0{\/N'^) structure of the /3- function here we are putting in 
place a means of checking any future attempt at that difficult programme. 
However, more importantly, our motivation in the present calculation goes 
beyond the two dimensional model for two reasons. First, we are determin- 
ing the d-dimensional structure of the /3-function by exploiting the conformal 
nature of the fixed point and applying the d-dimensional generalization of 
the star-triangle relation, we are able to deduce new information on the three 
dimensional model simultaneously. Our result goes beyond anything which 
has been calculated before in this area since we are computing at 0{1/N'^) 
and this is relevant for recent lattice simulations of the three dimensional 
four fermi interactions for relatively low values of N which have been a topic 
of interest lately, [^] . 

Further, the anomalous dimensions of all the fields of the theory have 
been deduced in 1 1£ ] and by calculating A we will in fact solve the model com- 
pletely to second order since the remaining critical exponents characterizing 
the model can be deduced through the hyperscaling relations which have 
been checked at 0{1/N) in |2C] although their consistency merely reflects 
the renormalizability of the model. Indeed the work we present builds very 
much on the calculation of A in the 0{N) case, Although there are sub- 
stantially more Feynman graphs to consider due to the presence of another 
field it is possible to make use of results obtained for those graphs in the 
continuous chirally symmetric case to minimize the work needed to obtain 
new results. This is significant when one realises the basic intractability that 
a conventional large analysis would entail. The previous remarks are im- 
portant in stating our second motivation on the need to have d-dimensional 
results. One feature of four fermi theories with various chiral symmetries is 
that they are believed to lie in the same universality class as the Yukawa 
models with the same chiral structure. Such equivalences, for example, are 
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established using the technique of the e-expansion. In essence the renormal- 
ization group functions, calculated to several orders in perturbation theory, 
are expanded in powers of e, where d = 2 + eor4 — 2e depending upon 
which dimension the coupling constant of the model is dimensionless in, and 
compared with the analogous function in a similar expansion in the other 
model. As the critical exponents calculated within the large N self consis- 
tency approach discussed here encode all orders information on perturbation 
theory, such d-dimensional results have been important in confirming such 



equivalences in the 0{N) case, |18 |. 



The paper is organized as follows. In section 2, we review the basic 
formalism we need for solving the chiral Gross Neveu model and recall the 
relevant features of the 0{\/N) calculation of A, which unlike the bosonic 
0{N) a model calculation at the same order requires the inclusion of several 
2-loop Feynman graphs. We devote section 3 to the derivation of the master 
equation whose solution will yield A at 0{1/N'^). This involves considering 
several three and four loop graphs in the Dyson equations of two of the fields 
of the model. The equation is solved in section 4 to obtain an arbitrary 
dimensional expression for the exponent from which we deduce information 
on the 3-dimensional model simple corollary. 

2 Preliminaries. 

The lagrangian of the model we will analyse is, 

L = zVS'#' + (TV5V' + «7rV^S^^'-^(f^^ + ^^) (2.1) 

2c/^ 

where a and tt are auxiliary bosonic fields which when eliminated from (2.1) 
yields a model with a four fermi interaction which has a [/(I) x U{1) global 
chiral symmetry. The fermion field ip^ lies in an A^-tuplet, 1 < i < N, 
and the parameter 1/A^ will play the role of the coupling constant in our 
analysis and we use the convention that tr 1 = 2. In the usual approach 
to solving (2.1) in the large expansion, [|l|, one can perform the path 
integral over the fermion fields in (2.1) as they appear quadratically in the 
action. To proceed further one uses the saddle point approximation when N 
is large to deduce from the effective potential that the true vacuum of the 
theory is not the perturbative one but that where the a and vr fields become 
dynamical. Further, in this vacuum a mass is generated for the fermion 
field although perturbatively it remains massless, ||l|. Equipped with this 
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structure one normally proceeds to solve the model in the large expan- 
sion by renormalizing the various Green's functions. Consequently, one can 
deduce the large approximation of the /3-function, for example. To go 
beyond the leading order turns out to be impossible due to the appear- 
ance of intractable integrals whose divergence structure, handled by some 
regularization, cannot be extracted. This is chiefly a consequence of the 
mass present in the propagators of the ip and a fields where the form of the 
latter involves a hypergeometric function of the momentum and mass, 
(This non-fundamental structure is consistent with its fermion bound state 
nature.) 

As we have already mentioned in the self consistency approach of ||2| it is 
indeed possible to go beyond the leading order by examining the structure 
of the field theory in the neighbourhood of its d-dimensional critical point 
defined as the non-trivial zero of the /^-function. There the theory possesses 
a conformal symmetry and the Green's functions of the fields obey a sim- 
ple power law behaviour and there is no mass in the problem. By solving 
for the critical exponents of the power law one can deduce information on 
the renormalization group functions through a critical point analysis of the 
renormalization group equation which takes a solvable and simplified form at 
criticality, For this paper, taking as the coupling constant allows us 
to compute the exponent order by order in by solving the appropriate 
Dyson equations truncated to the order we are interested in. 

We will now review the 0{\/N) derivation of 2A = — (3'{gc) in the self 
consistency formalism to illustrate the general method and to define the 
notation of the problem. As the a and vr fields are dynamical, we take. 



B 



^(^) - 7^[l + (^')'^'] (2-3) 



nix) ~ -£-[! + (x')^C'] (2.4) 



as the asymptotic scaling forms of the propagators of the fields in coordinate 
space in the critical region as x — > 0. The quantities A, B and C and A' , 
B' and C are the respective x-independent amplitudes, whilst a, f3 and 7 
are the full dimensions of the respective fields. We choose to compute in 
coordinate space in order to implement various integration techniques which 
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are easier to apply there but note that one can access the momentum space 
formulation easily via the Fourier transform 



valid for all a where a(a) = r(/i — a)/T[a) and we set the spacetime di- 
mension d to be d = 2/i for convenience. One can deduce the canonical 
dimensions of a, [3 and 7 from a dimensional analysis of the action but 
quantum fluctuations through, say, radiative corrections will always adjust 
these. To allow for this one introduces the anomalous dimensions in the 
definitions of the exponents via 

a = ^x + \ri , /3 = l-rj-Xa , 7 = ^ - V - Xtt (2.6) 

where rj is the fermion anomalous dimension and Xo- and Xtt are the anoma- 
lous dimensions of the 3-vertices of (2.1) 

In [0, |l^ we considered only the leading order terms of the asymptotic 
scaling forms in the skeleton Dyson equations which allows one to deduce 
the exponent t] at 0{1/N'^). To determine /3-function corrections one has to 
include the higher order corrections as we have in (2.2)-(2.4) which is similar 
to previous analyses in other models, |^ 0]. Then to deduce the exponents 
one writes down a set of equations which represent the Dyson equations in 
the critical region and solves them self consistently. We have illustrated the 
Dyson equations we need to consider for Ai, where X = fi — 1 + J2i^i ^i/^'^-, 
in fig. 1. Ordinarily to compute, say, 771 one truncates the equations at one 



loop, However, as was noted in |21] for field theories where the basic field 



is fermionic one has to include the two loop graphs of the a and vr 2-point 
functions where each has an (x^)^ insertion from (2.3)-(2.4) on the cr or vr 
internal propagator. The reason for their inclusion is due to a reordering 
of the powers of N in the master equation whose solution yields Ai. The 
quantities V'"^! and vr"^ in fig. 1 are the 2-point functions and their 
asymptotic scaling forms can be deduced from (2.2)-(2.4) by inverting each 
expression in momentum space via (2.5) and they take the x-space forms. 



^-\-) - :4^§)l^[l-^'-(« (2-7) 
^-Hx) - (2.8) 
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^-'(^) - cW^i^-C'Q{l){x')'] (2.9) 
where we define 

a{a — ^) a(a — /U + \)a{a — A) 

pya) = — 5 — -—^ , qia) = 

7r^'^a(a) a(a — iJ,)a[a) 

(/i — Qj (a — Aj(a — ;U + Aj 

If we formally denote the values of the 2-loop integrals by IIia,b,c and 
^ia,b,Cj where the subscript A, B or C denotes an insertion on a V', c or vr 
line respectively, then with (2.2)-(2.4) the Dyson equations of fig. 1 become 

= r{a-l)[l-A's{a-l){x^)^] + z[l + {A' + B'){x^)^] 

+ y[l + {A' + C'){x^)^] (2.11) 

^^^^ [1 - B'q{p){x'^)^] + 2z[l + 2A'{x^)^] (2.12) 



N 

- z^[Ui + {x^)\A'UiA + B'Uib)] + zy[U2 + ix^)\A'U2A + CUc)] 
= ^[l-CV7)(x')'] + 2y[l + 2A'(x2)^] (2.13) 

- y^[Ai + {x^)\A'A,A + C'Aic)] + zy[A2 + {x^)\A'A2a + B'A2b)] 

where we have set z = A^B and y = A^C. Comparing powers of x'^ the 
equations decouple into a set which become the consistency equations for r]i 
and give, |l|], 

2r(2/i - 1) 

r(^-i)r(i-M)r(^)r(^ + i) ^ ' > 

and a set involving A' , B' and C . The solution of the latter yield Ai where 
the relevant consistency equation is formed by setting the determinant of 
the matrix formed by treating A' , B' and C as independent basis vectors, 
to zero. In this matrix one can examine the dependence of each element 
and observe that r{a - l)s{a - 1) = 0(1) whilst p{l3)q{(3) = 0{l/N). Since 
y and z are both 0{\/N) from the r/i equation it is easy to see the necessity 
of including the higher order graphs IiiB,c and Aib,c- Direct calculation, 
however, reveals that Hib = Hic = Aib = Aic = 11, where the basic integral, 
n, has been computed explicitly in the 0{N) model, 21]. It turns out 
that in manipulating the determinant through elementary row and column 
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transformations and using, |19] 



- « - ^ (^-15) 

that one obtains 

' s{a-l) 1 1 

det I 2 q{f3) zU j (2.16) 

g(7) - zn 

where the Il-type contribution has cancehed in the central element through 
the additional symmetry present, compared with |l9|. Consequently, the 
consistency equation which results from (2.16), 

= q{P)s{a - 1) - 2 (2.17) 

does not require the value of 11 and we deduce that 

Ai = - (2^ - l)r/i (2.18) 



which is in agreement with [20| and establishes the correctness of our proce- 
dure for this model. This is a non-trivial statement as a second consistency 
equation would appear to emerge from (2.16) by setting the lower right el- 
ement to zero. We discard this as it does not give a result consistent with 
the exponent calculated in the conventional large N approach. Moreover, 
the appearance of such additional solutions is an artefact of the formalism 
and has been observed in other models where again it was appropriate to 



disregard it, [22|. Indeed from the transformations of the determinant we 
have made this potentially alternative solution does not take into account 
any contribution from the ip consistency equation, which again justifies our 
stance. 



3 Master equation. 

To proceed beyond (2.18) and to compute the 0{1/N'^) correction requires 
that one expands the quantities s{a — 1), g(/3) and 5(7) to the subsequent 
order and also includes the higher order corrections to the Dyson equations. 
For the former this requires the values of the vertex anomalous dimensions 
and 7/2- They have been calculated in [19| by considering the subsequent 
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corrections to the Dyson equations which determine the fermion anomalous 
dimension as 



V2 = Vi 



(3.1) 

X.i = Xni = (3.2) 



where ^'(/u) = ip{2fi—l) — ^(1) + Tp{2—fi) — ')p{fi) and ^{x) is the logarithmic 
derivative of the T-function. By considering the scaling behaviour of the 3- 
vertices at 0{1/N'^) one also finds that 

V - V - ^^(V-10/i + 7)r?^ 

X.2 - X.2 2(^-1)3 ^^-^^ 

using methods developed from the earlier work of p3| , p^ j. 

The main problem in determining A2, however, is the inclusion of the 
higher order corrections. For the fermion equations this is straightforward 
in that one includes the two loop Feynman graphs of fig. 2, and these 
together with the two loop graphs of fig. 1 are all that is required to de- 
duce r]2 from the leading order ansatze of (2.2)-(2.4) and (2.7)-(2.9). It is 
worthwhile detailing the resulting Dyson equation for ■0 explicitly since it 
will illustrate several important points which arise in the a and tt equations 
and may be obscured there. We have 

= r{oi-l)[l- A's{a-l){x^f] + z{x^Y^[l + {A' + B'){x^f] 
+ y{x^Y-[l + {^ + C'){x^f] 

+ z2(^2)2x.p^ + (x2)^(A'SiA + 5'Sib)] 

- 2yz(x^Y^^^^ [S2 + {x^)\A'J:2A + B'^2B + C"S2c)] 

+ y\x^f^-[^^ + {x^)^{^T.^A + C'^zc)] (3.4) 

where the mixed two loop graphs have the same value and we have not 
cancelled the powers of x^. The values S and Sa,_b,c correspond to the 
values of the integrals representing the Feynman graphs without symmetry 
factors which have been displayed explicitly. In analysing similar equations 
at 0{1/N'^) in other models the analogous higher order two loop graphs are 
infinite which can be deduced by a detailed computation. To handle such 
divergences one shifts the exponents of the a and vr fields by an infinitesimal 
amount A through Xo- ^ Xo- + ^ ^iid x-w ^ X-w + ^ where A plays the role 
of a regularizing parameter. The subsequent simple poles in A which would 
occur in the regularized graphs are renormalized by the vertex counterterm 
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available from the one loop graphs. For the current equation it turns out 
that upon adding the divergent graphs together the infinity cancels due to 
the symmetry induced through the 7^ interaction and there is therefore no 
need to renormalize (3.4) to the order we are working to, 0{1/N'^). Thus 
only the A-finite parts of the 2-loop graphs remain in (3.4). As before we 
decouple the finite equation into that which has already been used to deduce 
r/2 and 

= [z + y-r(a-l)s(a-l) + z2SiA-2yzS2A + y^S3A]^' 

+ [z + z2SiB-2yzS2ij]S'+[y + y2s3c-2yzS3dC" (3.5) 

where the powers of are absent at this order due to (3.2). The explicit 
evaluation of the finite parts of Sj^ and Sjc*, however, have been determined 
in Q and it turns out that to the order we are working to they are zero. 
Also if one compares the iV-dependence of each term of the coefficient of A' 
in (3.5) it is easy to see that T,iA is 0{1/N'^) relative to r(a — l)s{a — 1) and 
therefore can be neglected in the correction to the determinant of (2.16). 
This term would be relevant only for A3. Thus the part of (3.5) which we 
require for A2 is simply 

= [z + y-r{a-l)s{a-l)]A' + zB' + yC' (3.6) 

where, of course, the 0(1/A^^) parts of z and y contribute in the last two 
terms. 

At leading order we had to include the higher order two loop graphs of 



fig. 1. As was noted in [21| this is a feature of the formalism for models 
where the fundamental field is fermionic, although in the final equation, 
(2.17), for the chiral Gross Neveu model the explicit value was not needed. 
By the same argument one has now to include the representative three and 
four loop higher order graphs of fig. 3 in the a and vr equations which we 
now discuss. We have given only the basic distinct topological structures of 
the graphs which occur for a where the internal bosonic lines can be either 
the (T or vr fields. So, for example, there are 4, 4, 4, 8 and 8 respective 
graphs for each of the basic ones given in fig. 3 for the a equation. Of 
course, since tr7^ =0 only half of the four loop graphs survive. Further, 
when one substitutes the asymptotic scaling forms (2.3) and (2.4) one need 
only include those graphs in the a consistency equation where there is one 
(x^)^ insertion on an internal bosonic line, in a similar way to leading order, 
P]. The full consistency equation for a with the graphs of fig. 3 is quite a 
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long expression and rather than give its complete form we write down only 
the part relevant for A2 as 



= 4zA'-B' 
where we have set 

nc2 = 



N 



+ C'lyzU - Uc2] 



(3.7) 



2(n2iji + n2B2) + n3B +n4s 

2(2n5Bi + U5B2) - m^Bl + 2UeB2) 
2{U2C1 + n2C2) + + II4C 

2(2n5Cl + n5C2) - 4(n6C71 + 2n6C2) 



(3.8) 
(3.9) 



and our notation is partially defined in fig. 3. When there is an additional 
subscript 1 or 2 in (3.8) and (3.9) this corresponds to the two distinct ways 
of including an insertion on the internal bosonic line and is consistent with 
definitions given in Using the properties of the 7-matrices and being 
careful in manipulating factors of i from the vertex of (2.1) one finds that 
the contributions from 1X251, ^3B, ^2Ci and Use style graphs each sum to 
zero through use of (2.15). This is consistent with the fact that one does 
not require a regularization at this order in 1/A^ as the divergent contribu- 
tions cancel. Subsequently this implies Xai = which is in agreement with 
calculations of other Dyson equations. Further, with (2.15) it is possible to 
simplify (3.8) and (3.9) to 



Ii-B2 

nc2 



4z'[U2b2+K 



ttO I ttO 

~''^2B2 + -"-MB 



zi(2n°Bi+n°52 
-zi(2n°5i + n052 



+ 2n°Bi)] 



&Blj 
QBl) 



(3.10) 
(3.11) 



using (2.15) where the superscript ° denotes the value of the fundamental 
graph, with the appropriate insertion, given in fig. 3. Whilst these have all 
been computed explicitly in the 0{N) case, we defer to later the direct 
substitution into 1152 and Ilc2- 

A similar analysis of the vr Dyson equation yields the third of our con- 
sistency equations as 



= Ay A' + B'[yzU - Uc2] - C 



P(7)g(7) , 2rr , TT 
j;^ h y n + IIb2 



(3.12) 



where the same combinations 1152 and Ilc2 appear as in (3.7). We can 
now manipulate these formal corrections to (2.11)-(2.13) to deduce the mas- 
ter equation for A2. First, though we record that from the 772 consistency 
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equation, |1£], 



Z2 



(3.13) 



and following the same transformations of the determinant as we made at 
leading order with these corrections we are left with the determinant of the 
2x2 submatrix 







det 



2z 



r{a — \)s{a 
Az 
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2z 



(3.14) 



Again the two loop graphs have cancelled though it is important to note that 
they would in fact have contained 0{1/N'^) information since the integral 
depends on the exponents ryi and Ai and they would have contributed to A2. 
From (3.14) we deduce that A2 will emerge from the solution of 



[r{a - l)s{a - 1) - 2z\ 



N 



+ IiB2 + 



(3.15) 



Again we ignore the solution that would come from the lower right element 
of the full 3x3 for the reasons stated earlier. Moreover, we note that the 
higher order graph contributions appear in a particular combination 11^2 + 
11(72 and from (3.10) and (3.11) this simplifies to 



nB2 + nc2 = 8z3[nOB-zi(2n°5i + nOB2)] 



(3.16) 



Thus we need only consider the contributions from three of the eight basic 
topologies and insertions of fig. 3 which represents a significant simplifica- 



tion. We note that the graph n^^^ corresponds to an 
top or bottom internal bosonic line. 



„2\A 



insertion on the 



4 Discussion. 

All that remains now is the explicit evaluation of (3.15) at 0(1/A^^). The 
values of the higher order graphs have been computed in Q for the Gross 
Neveu model with a discrete chiral symmetry. This involved extensive use 
of the technique known as uniqueness, |p, which is applicable to calculating 
massless Feynman diagrams, and is a method of integration which is based 
on the conformal properties of the integral in d-dimensions. Our manipula- 
tions for (2.1) have been such that the basic graphs of (3.16) are equivalent 
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to the values of the graphs of the 0{N) case and therefore we record that, 



vr 



4fM 



(A.-i)2r4(^) 



39 + 



(4.1) 



5B1 



7r6A'(2/i - 3)a(2^ - 2) 
2(^-l)5(^-2)r3(;x) 



69 _ $ - ^r^ + 



2(/.-l)2 



2(u - 2)^- 

+ ^7^^ TT- + 



(2^-3)(^-i) ' 

{fi-2) 



(2/i-3)(M-l) 



+ 



(2^ - 3) ii^- 1)2 



(4.2) 



^hB2 



(/i-2) 



+ 



7r^'^a(2/i - 2) 

(3;U-4)^ 
(^-l)(/.-2)2 ' (;,-2)2 



+ 



2(/i-l)V 
27r6Ma(2/x - 2) 
(/i- 1)6(^-2)2 



(4.3) 



where G(/i) = V'(/x) - V'(l) and ^>(//) = V'(2/i - 1) - ^'{2 - /i) - V''(m) + 
■(/''(I). Thus substituting into (3.16) we have 



IIb2 + 



8(^-1)^ 
/x(/i - 2)2r/i 
2(2^-3) 



2(2^ - 3) 



3(3^-4)9 
(/^-2) 



2(2;U-3)(/i2 -4^ + 5)^- 



+ 



(A.-l)(^-2)2 

6 1 

+ 



+ 



(4.4) 



(^-2)(^-l)2 (;,-2) (^-2)2. 
A httle bit of algebra subsequently leads us to the main result of the paper, 
3;u2(3;U-4)9 



A2 



;u2(2^-3)(cI> + M/2) ^ 4/x(/i-l) 



2(/x-l)(^-2) (^-l)(;,-2) 



if, - 2)2r?i 



+ 



+ 



1 /i2(2^ - 3)(2/i3 - 9/^2 + 16^ - 11) 



(/i-1) ■ 2fx 
+ 2/i(2/x - 3) + ^' 



2(^-1)3(^-2)2 
V2(2^-3)(;u2 -4/i + 5) 

(/X- 1)2(^-2)2 



4/i2 + 1 



(4.5) 



in arbitrary dimensions. 
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As the three dimensional chiral Gross Neveu model has also been the 
subject of study recently, [18, 2C, we can set /x = | in (4.5) to deduce a 
new result for this case as 

1 16 16[472 + 277r^] 
^ - 2 " 3^ + ThT^ ^^-^^ 

This will be important for deducing numerical estimates for this exponent 
for relatively small values of N since several orders are now known. For 
example, in Q an estimate was given for the 0(8) Gross Neveu model which 
agreed with a recent Monte Carlo result, |l8|| . 

We conclude with various remarks. First, the chiral Gross Neveu model 
has now been solved completely at 0{1/N'^) and it is important to recognise 
that this is one order further than had previously been possible. By this 
we mean that two independent critical exponents A and r] ox r] + Xa have 
been determined to this order. The remaining thermodynamic exponents 
can be deduced through the hyperscaling laws which are valid for renor- 
malizable quantum field theories and which have been checked in |^0|. The 
quantum structure of the model is appreciably simpler than for the 0{N) 
Gross Neveu model in the sense that fewer topologies of Feynman diagrams 
needed to be considered and this is primarily due to the extra chiral sym- 
metry present in (2.1). Indeed a clear indication of symmetry within the 
self consistency formalism is the vanishing of the anomalous dimensions of 
one of the fields or vertices, such as Xo-i = XttI = here, in all dimen- 
sions. A similar phenomena occurs in the supersymmetric 0{N) a model 
where the anomalous dimension of the analogous supermultiplet of auxiliary 
and Lagrange multiplier fields vanishes at 0(1/N) as a consequence of the 
unbroken supersymmetry present in the model, pC 
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Figure Captions. 

Fig. 1. Leading order Dyson equations to determine Ai. 

Fig. 2. Higher order graphs for ip Dyson equation. 

Fig. 3. Higher order graphs for a and it Dyson equations. 
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